Solving for tadpole coefficients in one-loop amplitudes 
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Abstract 

One-loop amplitudes may be expanded in a basis of scalar integrals multiplied by rational coefficients. We 
relate the coefficient of the one-point integral to the coefhcients of higher-point integrals, by considering the 
effects of introducing an additional, unphysical propagator, subject to certain conditions. 



1. Introduction 

One-loop scattering amplitudes may be expanded in a sum of scalar integrals 0, [3, 0| multiplied by 
rational coefficients. This expansion arises explicitly in typical computational approaches, reviewed for 
example in 0] . The coefficients may be derived directly by reduction of Feynman integrals [H, 0] , or they 
may be sought as solutions to linear equations taken from various singular limits, such as unitarity cuts, in an 
on-shell formalism Q. Within the second approach, the coefficients can be found by applying "generalized 
unitarity" multi-cuts @, 0> S]- Alternatively, since the master integrals are known explicitly and feature 
unique (poly)logarithms, they can also be distinguished by the usual unitarity cuts, which are double-cuts 
0. One way to do this is by rewriting the measure of the cut integral in spinor variables, and then applying 
the residue theorem [lo| . 

This procedure of spinor integration has been carried out in generality for theories with at most scalar 
masses, and analytic expressions for the coefficients of the scalar pentagon/box, triangle, and bubble integrals 
have been given [lT| . However, the tadpole coefficients are missing, simply because they are obviously free 
of cuts in physical channels. Our note addresses this point. 

We find that we can solve for the tadpole coefficients in terms of the coefficients of higher-point integrals 
after introducing an auxiliary, unphysical propagator. The auxiliary loop integral then has two propagators, 
so we can apply unitarity cuts formally. The tadpole coefficient is accordingly related to the bubble coefficient 
of the auxiliary integral. Our result is a set of relations giving the tadpole coefficients in terms of the bubble 
coefficients of both the original and auxiliary integrals, and the triangle coefficients of the auxiliary integrals. 
It is interesting to consider whether this construction might have other applications. 

To derive the relations between tadpole coefficients and the others, we make use of work of Ossola, 
Papadopoulos, and Pittau (OPP) Jj, which gives the result of one-loop reduction at the integrand level. 



building upon analysis of their tensor structure In addition to the integrands for scalar boxes, triangles, 
bubbles, and tadpoles, there are a number of "spurious terms" which vanish after integration. The complete 
decomposition and classification given by OPP allows us to relate the original loop integral to the auxiliary 
integral including the unphysical propagator. We then derive relations among their respective coefficients, 
and identify conditions that almost completely decouple the effect of the unphysical propagator. 

We note that on-shell approaches to loop amplitudes face important subtleties in seeking tadpole coeffi- 
cients analytically. The operation of making a single cut relates an n-point loop amplitude to an (7i4-2)-point 
tree level quantity, which should be considered as an off-shell current. These are the same starting points as 



in proposals to reconstruct full amplitudes entirely from single cuts [13|, |14| . Another cut-free integral is the 



0-mass scalar bubble. In this note, we assume that this contribution is known. For applications to physical 
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amplitudes using unitarity methods, it will be necessary to account for cuts of self-energy diagrams [8|. The 
procedure for finding massless bubble coefficients is not known in general, analytically, and merits further 
exploration; here, we will comment on the possibility of taking a limit of vanishing mass. Another proposal 
[15i] is to fix the tadpole and massless bubble contributions by universal divergent behavior, once all other 
integral coefficients are known. 



2. Relations among cuts and coefficients 

We adopt the notation of OPP 0, ■ The D-dimensional loop momentum is denoted by g, whose 4- 
dimensional component is q. The denominator factors take the form Di — (q+pi)'^ — Mf , where i = 0, 1, 2, . . .. 
The tadpole of interest shall be associated to the factor with i — 0. We define £ = —q — pQ. and Ki = po —pi- 
Expanding the loop momentum variable into its four-dimensional component plus the remaining part q 
satisfying q^ — — the denominators can be rewritten as Di — {£ — K.if' — Mf — p? . 

We are interested in the effect of including an auxiliary denominator factor, which we write as 

Dk = {I- Kf - Ml^ - fi^. "'-p-p' (1) 

At this point, K and are variables unrelated to the physical amplitude. Later, they will be chosen 
subject to conditions that minimize the effect of this auxiliary factor. 
The one-loop integrand is 

N{q) 

[truc-intcgrand] / n\ 

UqUi ■ ■ ■ Um-1 

where, following 0], we use N{q) to denote the numerator, which is a polynomial in q. We call the integrand 
Itrue the "true" integrand to distinguish it from the "auxiliary" integrand, which we construct by inserting 
the auxiliary factor Dk, as follows. 

J- ^ (.Q^ [aux-integrand] /q'j 

~ DkDqD^ ■ ■ ■ Dm-l ^ ' 

Consider the single-propagator cut of the tadpole of interest. It is the result of eliminating the denomi- 
nator factor Z?o from the integrand: 

4 tree 

N{q) 

[l-cut-amp] 

Ul - ■ ■ Um-l 

This integrand is the analog of the product of tree amplitudes ^Lrf't^afght obtained from a standard unitarity 
cut. However, from a single cut, we obtain a tree amplitude at a singular point in phase space, since two 
external on-shell momenta are equal and opposite. This singularity can create difficulties that we do not 
address generally here. It is probably best considered as an off-shell current. 

In the OPP method j?!, the integrand is expanded in terms of the master integrals multiplied by their 
coefficients in the amplitude, plus additional "spurious terms" which vanish upon integration. The uninte- 
grated expansion is 

ru — 1 m— 1 m— 1 

Itrue - 5]Kz)+5((7;z)]/(')+ + ^ [c(z, j, r) + z, j, r)]/(^-^''-) 

i i<j i<j<r 

m — 1 m — 1 

i<j<r<s i<j <r<s<t 

where a(i), j), c(i, j, r), d{i, j, r, s), e{i, j, r, s, t) are the coefficients of the master integrals; a(g; i), h{q\ i,j), 
c{q; i, j, r), d{q; i, j, r, s) are the spurious terms which integrate to zero; and the master integrals are 

/(«) — J_ /(ij) — _ ^_ j(i,j.r) _ _ } _ j{i,j,r,s) _ _ _ j{hj,r,s,t) ^ 



Di D,Dj D^DjDr D,DjDrDs D.DjDrD^Dt 
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Notice that we have included the pentagon exphcitly. We shall perform our analysis of the coefficients in 
D — A — 2e dimensions. 

Now, consider the auxiliary integrand Ik- On one hand, it is simply the true integrand divided by the 
auxiliary propagator. Thus, using ([5]), we get the following expansion in master integrals: 

J m—l m — 1 m — 1 

i i<j i<j<r 

m—l m — l 

i<j<r<s i<j<r<s<t 

Notice here that the "spurious" terms such as / b{q; are no longer spurious with the factor Dk included. 
For example, while / 6(g; = by construction, in general f b{q; i, j)/'^^'*'^^ =^ 0. On the other hand, 

the auxiliary integrand Ik has its own OPP expansion, where we label the auxiliary coefficients and spurious 
terms by the subscript K, and we separate the auxiliary propagator explicitly, so it is not included in the 
summation indices: 

m — 1 m — l m — l 

Ik = Y^aK{i)+aK{q;i)]I^'^ + Y[^K{^,J) + bK{q;^,J)]I^'''^ + ^ [ck(*, J, r) + ^^^(q; j, r)]/'^'^-'') 

i i<j i<j<r 

m — l m — l 

+ Y [dK{i,J,r,s)+dK{q;i,J,r,s)]l(^'^^'''^^^+ ^ eK(*, J, r, s, t)/^^'^''-'^'*) 

i<j <r<s i<j <r<s,t 

m—l m — l 

+ [aK{K)+aK{q;K)]I^''^ + Y[bK{K,j) + bK{q;K,j)]I^^'^^ + Y{^K{K,j,r) +CK{q;K,j,r)]I^''^^'^^ 

3 j<r 
m—l m — l 

+ Y idK{K,J,r,s)+dK{q;K,j,r,s)]I^''^^'^^^^+ ^ eK(i^,J- (7) 

j<r<s j<r<s,t 

With the subscript K, terms such as 6a:(<z; hj) and bK{q', K,j) are truly spurious terms in e.g. / bK{q] K, j, r)l'^^'^^ 
0. 

The purpose of introducing the auxiliary integrand Ik is to give information about the tadpole coefficient 
by cutting two propagators. So, we will always choose to cut the auxiliary propagator Dk along with Dq. 
Restricted to the terms involved in this cut, the auxiliary integrand ([7]) is 

'm— 1 

Ik\c,, = [bK{K,Q) + bK{q-,KMl^''^''^ + Y.^ck{K,Q,i) + CK{q-,K,Q,i)]I^'''''^^^ 

i 

'm—l m—l 
i<j i<j<s 

Similarly, we can restrict our attention to the corresponding subset of terms in 

m — l m — l 

Ik\c,^ = [a(0)+a(g;0)]/(^^") + ^[6(0,z) + 6(g;0,*)]/(^'"'') + 5][c(0,i,j) + c^^^ 

i i<j 
m — l m — l 

+ Y ['^(O' J' ^) + ^(9; 0' r)]I^'''°^^'^^'^ + E ^(0' J' [a.v.ae._.„^^y, 
i<j<r i<j<r<s 

Our plan is to find the tadpole coefficient, a(0), by imposing the equivalence of ^ and ^ after com- 
pleting the cut integral. After integration, the spurious terms of (jS]) simply drop out, as they are designed 
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to do so: 



[ Ik = bK{K,0) [ /(^.") + X;'cK(i^,0,z) / 

'J Cqk Cqk j Cqk 

m—1 „ m — 1 „ 

i<j •^C'off l<j<S •''^OK 



Here, the cut integral is denoted by ^q^^i which indicates that we use the Lorentz-invariant phase space 
measure including the factor 5{Dq)5{Dk)- 

However, the integration of the formula ^ is not so straightforward, because the original spurious terms 
no longer correspond to the structures of the denominators they multiply. So, we shall view the expression 
Q as a function of the loop momentum q, and find the coefficients of master integrals analytically, for each 
of the spurious terms classified by OPP. 

Keeping in mind that our target is the single number a(0), which appears as part of the auxiliary bubble 
coefficient in ([3]), we begin by extracting only the auxiliary bubble contributions of the various spurious 
terms, divided by their denominators as well as Dk- (The other non-spurious terms, with 5(0, i), c(0,i,j), 
and d{0,i,j,r), clearly belong entirely to coefficients of other master integrals, of which the 4-dimcnsional 
pentagon is a linear combination of five boxes.) 

Our result is that there are conditions under which most of the spurious terms have no effect. Specifically, 
for all the propagator momenta Ki inside Di, we would like to take 

K ■Ki = 0, Vi; M|- = + . (11) 

We are free to take (fTT|) as a definition of while the condition for K is clearly nontrivial to satisfy 
physically. Let us consider this condition more concretely. For the purposes of defining our construction, we 
perform a formal reduction. Any integrand having five or more propagators has at least four independent 
momenta Ki that can be used to expand any external momentum vector appearing in the numerator to do 
the reduction. For integrands with at most four propagators, there are no more than three momenta Ki, 
and the condition pT|) can be satisfied, for example by the construction K^ — e^^p^K^Kl^K^ . 

If conditions PT|) are satisfied, then we find that only one of all the spurious terms contributes to the 
auxiliary bubble coefficient. Specifically, 

hK{K, 0) = a(0) + ^Y.^K^ - + Ml)bo^[Ki\, (12) 

i 

where &oo [Ki] is the coefficient of one of the spurious terms defined in Q . 

A convenient way to constrain &oo[-^i] is to identify the effect of the spurious term on the auxiliary 
triangle coefficient. Still imposing the conditions (fTT|) . we repeat our analysis of all the OPP spurious terms 
in ini, this time isolating the contributions to triangle coefficients. Fortunately, we find that only this same 
single spurious term has a nonvanishing effect, if we focus on the terms with /i^-dependence. (We assume 
that explicit /i^-dependence in the numerator N{q) has been set aside.) The result is 

ck{K,Q,i)\^, = h{Q,i)\^, + ^-boo[Kl (13) 

where |^2 means the coefficient of /i^. 

Now we propose the following procedure for finding tadpole coefficients. 

1. Find the single-cut expression A^l°1^^ obtained by cutting the propagator _Do- Expand the numerator 
in /i^, and work term by term, setting aside these explicit factors of . 

2. Construct the true integrand I — A*(^%^^/DQ and the auxihary integrand Ik = ^i-cut/(^if-^o)- It 
may be convenient at this stage already to choose K and Mk to satisfy the conditions pip . Alterna- 
tively, they can be taken as arbitrary variables until the final step. 

4 



3. Use the cut integral J(j^^ Ik to evaluate the auxiliary bubble coefficient bK{K, 0) and all the auxiliary 
triangle coefficients CK{K,0,i)- 

4. Use the cut integrals J(j^^ I to evaluate all the true bubble coefficients 6(0, i). 

5. The tadpole coefficient is given by imposing the conditions (jlip in the following expression. 

a(0) = bK{K, 0) + J2 ' 4j^2 ° i^M 0, ^) - 6(0, z)]|^. . (14) 

This formula is valid term by term, having set aside the original explicit factors of //^ in the numerator 
Niq). 



3. Contributions to bK{K,0) from spurious terms 



In this section we will discuss the contributions to 6if(i4r, 0) of the auxiliary integrand Ik from the 
expression where the terms have been separated into the scalar integral coefficients, plus spurious terms 
as classified by OPP fi\. As we have discussed, these terms are no longer "spurious" in the same sense, 
once Dk is included. (N.B.: OPP write the expansion with all denominators multiplied through, so that 
the "spurious terms" for them are the polynomial numerators. Here, we use "auxiliary spurious terms" to 
refer to the correesponding terms with all denominators present, including Dk-) 

The first contribution is obviously a(0), which is the tadpole coefficient that interests us. Now we discuss 
the possible contributions from the terms with a, b,c,d in © . We will see why we choose the decoupling 
conditions pT|) . They arise naturally by considering the terms of lowest degree. We have proceeded step 
by step through all the spurious terms of 0]. Our results have been derived in the formalism of [13] and 
verified using Passarino-Veltman reduction as implemented in FeynCalc [l6l |. 



• One-point spurious terms: In the simplest case, all spurious terms of this type are linear in the 
numerator In ~ 2£ ■ Ri/ Dq. The auxiliary integrand, including Dk in the denominator, is then 
/j^^ = 2£ ■ Ri/{DkDq). It is easy to find the scalar bubble coefficient from a standard unitarity cut 
(or alternatively, by straightforward reduction). The result is 

C[Dq,Dk\ = • (15) 

There are four independent "1-point like" spurious terms as given by OPP, i.e., four independent values 
of Ri . We see that we can decouple all their contributions by imposing the condition 

+ -Ml = ['"p'"'=»i (16) 

• Two-point spurious terms: Spurious 2-point terms can be either linear or quadratic in loop mo- 
mentum. In the case of linear dependence, the auxiliary integrand with Dk has the scalar bubble 
coefficient 

n 1 -{K ■ K^)iK ■ Ri) + KHK^ ■ Rl) [:™k1 

In the spurious terms, i?i takes three possible values of vectors, called l-j^ifi^n. These vectors are 
defined in here we only need to use some of their properties. (We use K as the auxiliary momentum 
in the OPP construction of these vectors.) In each of these cases, we have Ri ■ Ki = 0. Moreover, 
^■^7/8 — 0, but K -n ^ 0. To make this last spurious contribution vanish, we enforce a new decoupling 
condition: 
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(18) 



Now we move on to the quadratic spurious 2-point terms. There are five such terms. For four of them, 
the auxihary bubble coefficient vanishes under the two decouphng conditions. The fifth spurious term 
is K{q; 0, i), which can be written {i ■ n)^ — {{£ ■ Kif — Kfi'^)/3. Its coefficient in the OPP expansion 
is denoted 600 (0,i). After imposing the decouphng conditions, the auxihary bubble coefficient from 
this term is 

Because this spurious term gives a nonzero contribution under the decoupling conditions, we must 
calculate it and subtract its contribution when we calculate the tadpole coefficient. For this reason, 
we will turn to the auxiliary triangles ck^K, 0, i) in the following section. 

• Three-point and four-point spurious terms: All of the auxiliary three-point spurious terms 
decouple after imposing There is just one auxiliary four-point spurious term, and it gives no 

bubble contribution at all, because its numerator is linear in the loop momentum. 

To summarize, have seen that if we impose the conditions (jlip. then all contributions from spurious 
terms will decouple, except one, whose coefficient is 600 ■ We have 

bKiK,0) = aiO) + J2booiO,i) ' ^ ^° ^ (20) 

where &oo(0, i) is the coefficient of the spurious term K{q; 0, i) as defined by OPP. 

In this analysis, we are assuming a renormalizable theory. We have assumed that the power of £ in the 
numerator^is equal to or less than the number of propagators in the denominator. In those terms where the 
power of £ is strictly less than the number of propagators, then we have &oo(0, i) = 0, Thus we have 
bK{K,0) = a(0), i.e., we get the tadpole coefficient a(0) immediately by calculating the bubble coefficient 
under the decoupling conditions. From terms where the power of £ is equal to the number of propagators, 
&oo(0, i) 0, and we need to compute it. We have found that we can use a similar decoupling approach to 
calculate the triangle coefficient C[Dq, Dk , Di] and extract the corresponding 6oo(0,«). This procedure will 
be discussed in the next section. 



4. The calculation of boo(0,i) 

Recall the expansion ©, where we augmented the OPP expansion with the extra factor Dk, so that the 
spurious terms no longer integrate to zero. We see that the term b{q; 0, i)/(^'0'*) contributes not only to the 
coefficient of the bubble but also to the coefficient of the triangle /(^^O'*). Thus it is possible to find 

&oo(0,i) from the evaluation of the coefficient of triangle /(^^o^*) within Ik\cok- 

Just as in the previous section, where we studied all contributions to coefficient of 7^^'°^ from Ik\coki so 
CKiK,0,i) in ^ also receives contributions from the original spurious terms in Thus we carry out the 
corresponding analysis in this section. We impose the decoupling conditions from the start. Then we find 
that three of the auxiliary spurious terms still give auxiliary triangle contributions. The first nonvanishing 
contribution comes, indeed, from the term we want, namely K(q]0^i). Its auxiliary triangle contribution, 
after having applied the decoupling conditions, is 



The second and third nonvanishing contributions come from the spurious 3-point terms with quadratic 
dependence on the loop momentum. The auxiliary integrand is 



where i?i takes two values, called £3^4. After applying the decoupling conditions, we find that the triangle 
coefficient is 

{K ■ R,)mK, ■ K,){Kf + Mg - Mf ) ^ KfjKf + Mg - Mf)) 
K^{(K,-K,f-K^K]) 

This quantity does not vanish identically, so our decoupling might seem to be inadequate. But there is good 
news here: the contribution does not depend on fi^, while the contribution from 600 does depend on 
Thus we can use the /i^-dependence to find exactly the term we need. 

Our plan is now clear: (1) calculate the bubble coefficient &(0, i) from the integrand Itrue, which does not 
contain Dk\ (2) calculate the triangle coefficient CK{K^Q,i) = C[Do, Dk, Di] of the integrand Ik, which 
does contain Dk', (3) find the /i^-dependent terms in b{0,i) and Ci<-(A', 0, «), and then solve the following 
equation to find 600 (0, «). 

ck{K, 0,i)\^2 = biO,i)\^2 + ^boa{0,i), (23) 

where |^2 means the coefficient of ^i^. After computing 6oo(0, i) for every i, we substitute back into ([20]) and 
finally find a(0), the tadpole coefficient. 



5. Examples and discussion 

Finally, we list some formulas wc obtained from our algorithm and comment on their properties. We 
denote a general integrand term by two indices n, m, writing In,m[{Ki}, {RjY\ = YYj=i{'^^ ' ^j)/ (IliLo^ Di). 
Further, we define ai — Kf + A/q — Mf and A^j — {Ki ■ Kj)^ — KfK^. The tadpole coefficients from the 
first few of these integrands are: 

a[DoK^ = (24) 
a'[L>o\i22 = 3{K^)2 ^ ' 

rn 1 J:l=,A.AK.-R,)iK,-R2) 
a[Do]r,. = j^2]^^ . (26) 



where 



where 



Ai, ^ kI{Ki ■ K2), A22 = Kl{K^-K2), Au^ A2i = -KIkI l"^-*-^' (27) 



E,<fc A^^jsjK, ■ R2)(K, ■ R^){Ks ■ i?i) ,„3^^^, 



, rr-i 1 ^j^K -^J^y J - 133-tad I OQ\ 

A^0\l3i = olT^2 T^2\2\2^ ' (.28) 



= -2(^-1)' (Ai2(4qi(A-i • K2) - 2a2Kl) + 5K!{Ki ■ K2){a2{K^ ■ K2) - a^Kl)) 
^2;ii = A^,r2 = -2Kl{Klf{a^d.u-^Kl{a2{Ki-K2)-aiKl)), 

A2-22 ~ ^1;11 [ai ,ifi«a2,-K2 , j42;12 = ^1;22 = j4l ; 12 | ai , if i «Q,2 , ^2 



, Eti Et<.,i MtsjK, ■ R2)(Kt ■ Ri){K^ ■ fli) 
"t"""^^- = 'K!K!KiA,2A,sA2sA,2s ' ^''^ 
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where 

Ai23 = kIkIkI + 2(ifi ■ K2){K2 ■ KsKKs ■ Ki) - kI{K2 ■ Ksf ~ kI{K^ ■ Kif - Ki{Ki ■ Kif 

= 2/^23KIkI {(K2 ■ K3)[{K^ ■ K2fAis + {Ki ■ K-sfAi2] + {Ki ■ K2){K^ ■ Ki)\KlAri + i^a'Aia]} 

Ai;i2 = A2-ii=2K^^KlKiAi3A2-si~Kl{Ki-K3) + {K2-K^){K2-K3)), ^i;i3 = ^3,11 = [^i;i2]k2«K3, 

A2;22 = |ifi*-»X2 , — [Ai-ii]\ki^K3, j43;23 = ^2;33 = [j4i;12] |ki«K3 , ^2;12 = j4i;22 = [^1;12] |if i «/f2 , 

^2;23 = ^3;22 ~ [Ai;12]\ki ,K2,K3^K2,K3 ,Ki , As-^is — — [Ai-i2]ki ,K2 ,K3^K3,Ki ,K2 , 

Al-23 = ^2;13 = A3;12 = -2J^? 7^1 T^g^ Al2 A13 A23 

We note some patterns in these tadpole coefficients: (1) the tadpole coefficient is independent of /i^; (2) for 
I„^n-i, the coefficient is independent of masses; (3) for In,m the coefficient is of the form aiCi. 

Finally, we offer a comment on massless limits. In our calculation, we have assumed that Kf ^ 0. What 
happens when Kf 0? Let us use I21 as an example. The full reduction of I21 is 

Using the expressions for the scalar integrals, 

/2 (i^i = 0, Ml , M2 ) = (l-e)(Mg-M|) ' " ^ ^ ^ ° " ^ ' ^ ^ 

it is easy to see that the term with will be canceled exactly on the right hand side of l2,i- That is, for 

J2,i, knowing the coefRcient of the bubble, the smoothness of the limit Kf is enough to fix the tadpole 
coefEcient. Does this pattern hold in general? Looking back through our examples, we see that for J32, 133, 
and I43, the limits Ki,K2 — > give divergences which should be cancelled elsewhere. It is not yet clear 
exactly how generally we can relate the bubble coefficients to various tadpole coefficients by considering 
smoothness requirements in various limits, although the general behavior of Af ^ dependence suggests that 
such relations should exist. 
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